1

Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to “EPSG:3035” and check your changes.

You need to rename one of the id variables or adjust your join accordingly (“id = district_id”).

2

We want a first descriptive visual of the distribution of Covid-19 cases in Cologne and the surrounding districts. Calculate the number of Covid-19 cases in the last 7 days (cases_7days) by population (population) and multiply with 100,000 (in Germany usually called "7 Tages Inzidenzzahl).

Select Cologne (district_id == 5315), find the surrounding districts and plot the Covid Rate in Cologne and the surrounding districts.

You can use the dplyr function bind_rows to combine the two spatial objects “Cologne” and “Cologne Surroundings”.

3

Save your data set of Cologne and its surrounding districts as a ESRI Shapefile.